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1. Introduction 

Many nucleon structure observables require the calculation of disconnected quark line dia- 
grams for which all-to-all propagator techniques are needed. Here we present first results of an 
ongoing project to calculate the strangeness contribution to the spin of the nucleon \s as well as 
the scalar strangeness content of the nucleon (A/^l^^lA/^), using improved stochastic methods. 

The spin of the nucleon can be factorized into a quark spin contribution AZ, a quark angular 
momentum contribution Lq and a gluonic contribution (spin and angular momentum) AG: 

1 = 1ai:+l,+ag. (1.1) 

In the naive SU{6) quark model, AT = 1, with vanishing angular momentum and gluon contri- 
butions. In this case sea quark contributions will be absent too and therefore there will be no 
strangeness contribution As in the factorisation, 

AJ: = Ad + Au + As + --- , (1.2) 

where in our notation Aq contains both, the spin of the quarks q and of the antiquarks q. Ex- 
perimentally As is usually obtained by integrating the strangeness contribution to the spin structure 
function gi over momentum transfers x. The integral over the range in which data exists (x > 0.004) 
typically agrees with zero which means that a non-zero result relies on the unprobed very small-x 
region and is model dependent. Recent Hermes analysis [|l]] yields As = — 0.085(13)(8)(9) at a 
renormalization scale jj.^ = 5GeV^ in the MS scheme while our (as yet unrenormalized) results 
suggest lA^I < 0.01. 

The scalar strangeness density is not directly accessible in experiment but plays a role in 
models of nuclear structure. It is also of phenomenological interest since, assuming that heavy 
flavours are strongly suppressed, the dominant coupling of the Higgs particle to the nucleon will 
be accompanied by this scalar matrix element. 

We will first discuss our methods, then the error reduction achieved in our present lattice setup 
and finally we present results on the two matrix elements, before concluding. 

2. Stochastic methods 

We denote the lattice spacing by a and the lattice Dirac matrix by M = 1 — JC^. Disconnected 
quark line contributions require all-to-all propagators MJj^ where the multi-index / = {x, a, a) runs 
over all colours a = 1,2,3, spinor indices a = 1, ... ,4 and spacetime sites x (^V. Note that in our 
particular application it is natural and sufficient to restrict xto a given timeslice. Exact methods to 
obtain M^^ are unfeasible in terms of computer time and memory since 12V solver applications 
are required. Employing stochastic methods [^, this factor can be substituted by the number of 
estimates L <^ 12V: in a first step a set of Dirac noise vectors '■ i = 1, . . . ,L} is generated 
where the 12V complex colour-spinor-site components are filled with (Z2 /Z2)/\/2 uncorrelated 
random numbers These have the following properties: 

WjTIl ■■= TL\m){rii\ = 1 + , M = ^i^VI) . (2.1) 
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We will also employ the short-hand notation |-)(-| = |-)(-|^. We use the conjugate gradient algo- 
rithm with even/odd preconditioning to obtain the solutions \s^) of the sparse linear problems, 

M\si>) = \m>). (2.2) 

From these one can construct an unbiased estimate of ' : 

E{M-^) := \s){r]\ =M"^ +M"i (|t])(t]| - 1) . (2.3) 

V ' 

Due to the difference between E{M^'^) and M^^ above, any fermionic observable A can only be 
estimated up to a stochastic error Agtoch^ = ^( ^ / v^) on a given configuration. We define the 
configuration average over Wconf uncorrelated configurations and normalize this appropriately: 

OA,stoch •— — ■ U-^; 

"conf 

For large L and «conf this will scale like cr^stoch °^ (^"conf)^^- We also define the gauge error 
gauge °^ "conf Variation of the estimates of A over gauge configurations. This will be 

minimized at fixed ?iconf if A is calculated exactly. In general the gauge error is limited by, 

^A, gauge — *^/l,stoch • (2-5) 

If gjQ^j^ ~ a| gj^^jgg then obviously it is worthwhile to improve the quality of the estimates while 
stoch gauge ^hen precision can only be gained by increasing ?iconf > possibly reducing L to 

save computer time since the same n^^^^ scaling enters both sides of the inequality. 

In our calculation of A* the stochastic error initially was dominant. Hence we combined several 

variance reduction techniques to reduce this: 

• partitioning (also coined dilution) [Q]: we only set \r]i) / on one timeslice. This removes 
some of the (larger) off-diagonal noise elements, see eq. (E^, and reduces the variance. 



hopping parameter expansion (HPE) [^: the first few terms of the hopping parameter ex- 
pansion of Tr(rM^^) = Tr[r(l — JC^)^'] vanish identically but still contribute to the noise. 
For the Wilson action, Tr(rM"^) = Tr(rK-"^"M"') for « = 4, 8, depending on F, where for 
r = 1 one can easily calculate and correct for the zero-order difference. 

truncated solver method (TSM) |§|]: calculate approximate solutions \sn^/) after Ui solver iter- 
ations (before convergence), and estimate the difference stochastically to obtain an unbiased 
estimate of M~^: 



E(M ^) = K){r]\L, + {\s)-K)){n\L2 where L2«Li. 

Truncated eigenmode approach (TEA) |^, ||] : calculate the ?iev lowest eigenvalues and eigen- 
vectors of Q = 75M = Q\ Q"' = Qj' +L"fi \ui)qi^{ui\, and stochastically estimate the 
complement Qj' (with deflation included for free). 
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3. Lattice setup and error reduction 

Our exploratory calculations are performed onV = 16^ x 32 configurations of 2 + 1 
rooted stout-link improved staggered quarks with a Symanzik improved gauge action. These were 
provided by the Wuppertal group. The lattice spacing is fairly coarse, a^^ k, 1.55 GeV, and the 
spatial dimension is around 2 fm [^. We used the Wilson action for our valence quarks and 
currents with K = 0.166, 0.1675 and 0.1684, corresponding to pseudoscalar masses of about 600, 
450 and 300 MeV respectively. The analysis was performed on 326 configurations at Kjoop = 0.166, 
167 configurations at Kjoop = 0.1675 and 152 configurations at Kjoop = 0.1684, where K"ioop refers 
to the K value of the disconnected loop. Throughout we used a modified version of the Chroma 



code 

On each configuration the disconnected loop was calculated using the stochastic variance re- 
duction techniques detailed above (the TEA was only used at Kjoop = 0.1684, where 20 eigenvalues 
were calculated). We investigate the reduction in computer time, using optimized stochastic esti- 
mates, relative to those without any improvement techniques applied (except for time partitioning). 
We state all costs in terms of the average real computer time required on a Pentium 4 PC for one 
solver application (unimproved estimate), where we account for all overheads of the improvement 
methods. 



Tr(rioopM^i) 


'Qoop 


cost 


loop°P' 


^opt 

stoch 


loop 


''^stoch 


rioop = 5L;7i75 


0.166 


300 


-0.008(50) 


0.016 










100 


-0.033(55) 
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-0.185(148) 


0.135 






50 


-0.054(64) 
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0.1675 


300 


-0.085 (87) 


0.030 










100 


-0.040(101) 


0.054 
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0.198 
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-0.038(114) 


0.076 


0.056(265) 
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0.1684 


300 


-0.069(95) 


0.015 










100 


-0.068(96) 


0.036 


-0.089(216) 
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Tloop = 1 


0.166 


300 


14702.6(7) 


0.04 










12 


14702.5(7) 


0.18 


14703.5 (9) 


0.47 
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14702.3(8) 


0.23 


14703.7(1.0) 


0.65 




0.1675 


300 


14743.1(1.1) 


0.06 










12 


14743.4(1.2) 


0.33 


14745.0(1.3) 


0.69 






6 


14743.5(1.2) 


0.42 


14744.6(1.5) 


0.96 




0.1684 


300 


14764.9(1.2) 


0.04 










100 


14764.9(1.2) 


0.08 


14764.6(1.2) 


0.27 



Table 1: Results for the disconnected loop, averaged over configurations, obtained with (loop°P') and with- 
out (loop) variance reduction techniques. The cost is in units of the average computer time required to solve 
for one (undeflated) right hand side. 



Results for the configuration averages of the loops Tr(rioopM ^) are given in table |TJ The 
gauge errors CJgauge (that also depend on the stochastic noise) are displayed in brackets after the 



loop averages. These can be compared to the purely stochastic errors astoch> defined in eq. (2.4) 
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The deflation at K"ioop = 0.1684 where we apply TEA accelerates the solver but time is required 
for the eigenvector set-up. In our implementation the cost of solving for about 90 undeflated right 
hand sides equals that of 90 deflated ones (including this overhead). This is why in this case we do 
not display results obtained at the lower cost values. 

For Tr(^^y7y75M^') the stochastic error dominates over the gauge error unless L is chosen 
ridiculously large or variance reduction techniques are applied. Using these techniques the error 
is brought under control to the extent that we only need to invest the computer time equivalent of 
roughly 100 unimproved stochastic estimates to achieve astoch < ^Cgauge- In particular, we find a re- 
duction in CTg^toch (which is proportional to the amount of computer time required) of approximately 
25-30 for Kjoop = 0.166 and 0.1684. A smaller gain is obtained for the intermediate Kjoop = 0.1675 
which may benefit from using the TEA approach. For Tr(lM^^) the situation is reversed and the 
gauge error clearly dominates over the stochastic error: apart from possibly the heaviest JQoop there 
is no advantage in using variance reduction techniques. 

The matrix elements, 



{N,s\qy^yiq\N,s) = 2Mns^- 



(3.1) 



and {N\qq\N) are extracted from the ratios of three-point functions to two-point functions (at zero 
momentum): 



R^\t,tf) 



(r2pt<^2pt (^0 . ^/) Ix Tr(rioopM-' (x, f ; X, t))) 



(3.2) 



where r2pt = Tunpoi = (1 + 74)/2 and Hoop = 1 for (A^|^^|A^) and r2pt = /7;T5 (1 + 74)/2 and Tioop = 
7,75 for A^, where we average over 7 = 1,2,3. Note that for q = u,d there is an additional connected 
contribution 7?^^°", which we have not calculated. We combine the three Kjoop values with K"2pt = 
0.166 and 0.1675. In the limit of large times, tf^t^ to, 



{N,s\{qrioopqY^''\N,s 
2Mn 



(3.3) 
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Figure 1: The effective mass of the proton with K"2pt = 0.166 (left). The ratio, R^'>{t — 3a, tf) as a function 
of tf for fCioop = K-2pt = 0. 166 (right). 

We optimized the nucleon creation and annihilation operators using Wuppertal smearing with 
spatial APE-smeared parallel transporters [j8|]. The effective mass plot of figure [T] illustrates ground 
state dominance from a time t = 3a ^ 0.38 fm onwards. The same holds for K"2pt = 0.1675. Hence 
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we place the source at Iq = 0, the current insertion at ? = 3a and destroy the nucleon at > 4a. 
The result on the right of figure |l] does not depend on tf, even for tf < 6a, indicating that indeed 
with the chosen temporal separations we effectively realize the large-f limit. In table |^ we display 
the results for A^'*"* at the symmetric point tf = 6a ^ 0.76 fm: our methods enable us to reduce 
the squared errors by factors ranging from 5.5 to 11 at the fixed computational cost of 100 solver 
applications (in addition to calculating the two-point function). This falls somewhat short of the 
gains that we achieved in table |l] for the loops alone since now there are additional sources of gauge 
error. These we attempt to address in the near future. 





'Qoop ~ 
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+0.008(22) 


-0.004 (7) 


+0.008(20) 


50 


-^0.001(6) 


-^0.021(17) 


+0.004(10) 


+0.036(27) 






K-2pt = 0.1675 


300 


-0.005(6) 




-0.003(12) 




-0.004(13) 




100 


-0.008(7) 


+0.009(23) 


+0.005(15) 


+0.028(35) 


-0.006(13) 


-0.004(28) 


50 


-0.002(9) 


+0.046(29) 


+0.023(17) 


+0.083(51) 







Table 2: Results for Aq obtained with (R°P') and without (R) the use of variance reduction techniques. 



4. Results and Outlook 
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Figure 2: (A^|^^|A^)'^'^ (left) and A^'''^ (right) as functions of the quark mass used in the disconnected 
loop (expressed in terms of aMp^). The green points corresponds to a proton with K"2pt = 0.1675, while 
for the red points fC2pt = 0.166. 

In figure ^ we display our results for the two matrix elements where we obtained {N\qq\N)'^^^ 
at the cost of 12 solver applications per configuration and A^'^"* at the cost of 100 applications, in 
addition to the 12 applications that are necessary to calculate the two point functions. In neither 
case do we observe any significant dependence on the valence quark mass, varying this from nij^ w 
600 MeV down to 450 MeV, or on the loop quark mass, reducing nij^ ~ 600 MeV (~ strange quark 
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mass) to niji 300 MeV. We find |A5'| < 0.011 at the heavier proton mass and \^^s\ < 0.022 at 
the lighter mass value with 95 % confidence level while the scalar matrix element appears to be 
somewhat larger than one. Note however that the lattice results presented here are unrenormalized. 

In the near future we will further reduce the quark masses and the statistical errors, in partic- 
ular also of the scalar density, by refining our methods. We will also move to non-perturbatively 
improved Wilson sea quarks, allowing us to renormahze the results and to obtain a well-defined 
continuum limit. 
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